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Abstract 

We present computations of certain finite-size scaling functions and universal amplitude ratios in the 
large- limit of the CP^^^ field theory. We pay particular attention to the uniform susceptibility, the 
spin stiffness and the specific heat. Field theoretic arguments have shown that the long-wavelength de- 
scription of the phase transition between the Neel and valence bond solid states in square lattice S = 1/2 
anti-ferromagnets is expected to be the non-compact CP^ field theory. We provide a detailed comparison 
between our field theoretic calculations and quantum Monte Carlo data close to the Neel-VBS transition on 
a. S = 1/2 square-lattice model with competing four-spin interactions (the JQ model). 
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I. INTRODUCTION 

The emergence of relativistic gauge theoretic descriptions of complex condensed matter sys- 
tems at long wavelengths is an exciting theme of modem physics. It has now been established, 
as a proof of principle, that such descriptions can arise in suitably engineered lattice models of 
bosons'. The search for experimentally motivated microscopic models where such emergent phe- 
nomena may be observed is clearly of great interest, as this will pave the way to the detection 
and realization of these fascinating phenomena in experiment (e.g. in ultra-cold atoms in optical 
potential), and may hold the key to understanding the myriad mysterious properties of a growing 
number of strongly correlated materials. The most famous class of such materials are the cuprates, 
which begin as insulating 5 = 1/2 square lattice Neel ordered anti-ferromagnets, and then evolve 
into high-temperature superconductors on the introduction of a small density of charge carriers. It 
has long been held that a complete understanding of the insulating quantum anti-ferromagnet and 
particularly its transition to a paramagnetic phase from the Neel state holds the key to the cuprate 
mystery^. 

Motivated primarily by this notion, a large body of work on the paramagnetic phases of frus- 
trated, square-lattice, S = 1/2, SU{2) symmetric anti-ferromagnetic spin models has developed. 
Field theoretic work- established that the paramagnet that arises by condensing topological defects 
of the Neel order parameter is a valence-bond solid (VBS). Consistent with this field theoretic ar- 
gument, exact diagonalization of a S* = 1/2 square lattice anti-ferromagnetic spin model with a 
ring exchange on clusters with upto 40 spins have also found VBS phases on the destruction of 
Neel order^ as do series expansion studies on the Ji — J2^. There appears to be still some uncer- 
tainties on the Ji — J2 model which has been reviewed^. Recently, however an exciting step in 
this direction has been achieved: it has been possible to study the destruction of Neel order and 
appearance of the VBS phase on lattices with up to 10 000 spins*^ in an unbiased way in the so- 
called JQ model^ introduced by Sandvik by using sign-problem free quantum Monte Carlo (QMC) 
techniques: 



where the first term is summed on nearest neighbor bonds of the square lattice and the second term 
is summed over plaquettes allowing dot products only on nearest neighbor bonds. The JQ model 
harbors a Neel phase at J > 0, Q = and a VBS phase at J = 0, Q > and a transition between 
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them at J/Q ^ 0.040, Q > 0. In Ref.-^, the analysis of the QMC data provided evidence for a 
continuous transition. Subsequent QMC work has claimed a very weak first order transition using 
a "flowgram" analysis, although no detectable discontinuity in any physical quantity has been 
observedi^. Although the nature of this transition is currently under debate, further numerical 
work will likely be able to sort this out unambiguously, thanks to the absence of a sign problem in 
this model. Regardless of the fate of the transition at arbitrarily long length scales, there is clear 
evidenc o^i^i'° for very large correlation lengths and the associated scaling behavior on the relatively 
large intermediate length scales that have been simulated. In this paper, we take the natural point 
of view that this (possibly approximate) scaling behavior is the result of a nearby fixed point. An 
interesting challenge that then immediately arises is to identify the fixed point that gives rise to the 
observed scaling. 

A candidate for the fixed point was predicted in an extension^' of the field theoretic argument 
alluded to earlier. It has been shown that the long-wavelength description of the transition from 
the Neel state to the VBS state should be written in terms of the CP^ field theory of two complex 
bosonic spinons Za interacting with a gauge field, A^: 



with the constraint \za\'^ = 1, at = 2. An analysis of the Berry phases of the topological 
defects leads to the conclusion that only quadrupled monopoles of are permitted in the con- 
tinuum limit^''^. At = 1 duality transformations'^ establish that the quadrupled monopoled 
are irrelevant at the continuous transition of this field theory; they are hence also almost certainly 
irrelevant at the putative critical point at = 2. This leads to the remarkable conclusion that the 
long wavelength description of the Neel- VBS transition is described by Eq. ^ at N = 2 with 
a non-compact gauge field. It has been argued by Motrunich and Vishwanath^ that the critical 
point of the non-compact field theory at = 2 belongs to a new 'deconfined' universality class, 
distinct from the 0(3) universality class obtained when the gauge field is compact. Another study 
has also found evidence for the new universality classJ-^^. An accurate numerical estimate of the 
critical exponents and universal amplitudes of this new universality class are currently lacking, 
since Monte Carlo studies have been restricted to relatively small lattices. Numerical studies of 
both the non-compact^^ and compact'^ CP^^^ model are available. An alternate approach to a 
quantitative study, which we follow here, is to construct an expansion around the A^ = oo fixed 
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point. Indeed the N = oo fixed point is stable at large but finite A^, and universal quantities may 
be computed in a expansion and extrapolated to = 2. Our study assumes that there is 
a fixed point at = 2 that is continuously connected to the large- fixed point, as is true for 
instance in the 0{N) model; we however cannot prove that this is true. We note that our focus 
here is exclusively on the full SU (2) symmetric case; there has been extensive work^^ on a 
deformation of Eq. which was predicted^ to be the critical theory of a quantum transition be- 
tween a superfluid and a valence bond solid in U{1) symmetric spin models^, however our results 
do not apply to this case. 

The universality class of a fixed point is characterized by the values of the critical exponents, 
amplitudes and scaling functions close to the transition. We focus here on certain quantities that are 
associated with susceptibilities of conserved charges. We provide novel computations of (a) finite- 
size scaling functions for these quantities and (b) ratios of the universal amplitudes in the large-N 
expansion. We find reasonable agreement with QMC data on the JQ model, thus providing some 
support for the hypothesis that the scaling behavior observed close to the Neel-VBS transition in 
the JQ model is due to a proximity to the fixed point of the CP^ model. We note however that our 
calculations are only a first step and a fully convincing demonstration would require a comparison 
with a numerical study of an appropriate lattice discretization of the non-compact CP^ field theory, 
i.e. working directly at A^ = 2. 



II. LARGE-N FORMALISM 



We are interested in studying two-dimensional quantum anti-ferromagnets, that are described 
by Eq. at finite-temperatures. This clearly requires a study of the field theory in a slab geometry 
where the extent in the third direction is proportional to 1/T. We will be interested also in the 
effect of finite spatial extent of linear dimension, L. We describe the formalism used for these 
calculations in this section. 

We begin with the resolution of the constraint on Za, by introducing a real field A, which acts 
as a Lagrange multiplier at each point of space and time. 

^b = ^J d''rdT[\{d^-tA^)z^\^-tX{\z^\^-N)] (3) 

Note that the integration on r is carried out from to 1/T and the corresponding Fourier transform 
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consists of the Matsubara modes c<j„^ = 2TmrT, and likewise the spatial integral is from to L 
with corresponding Matsubara modes labeled by and Uy . In the limit N = oo the gauge field 
drops out and A takes on a uniform saddle point value. We can compute all quantities of interest 
from the free energy at criticality. We will organize its large-N expansion as, 

T = Nf + f^ + f^ (4) 



At = oo, the problem reduces to free complex scalar fields. 
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where we have included a convenient magnetic field H, which enables a computation of the uni- 
form susceptibility, Xu = 

We now turn to an overview of the computation of the free energy at next order, i.e. f^^ and 
/^^. We can organize the effective action for the A and A fluctuations as. 
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We will avoid explicit details of the computations of D12 and IIa here, since they have already 
been presented Ref.|20|. In terms of the functions D12 and 11^ the free energy is. 
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These expressions are useful to compute the corrections to the Wilson Ratio in Sec.[ 

The focus of this paper is on a computation of universal amplitudes associated with suscepti- 
bilities of conserved quantities. Before turning to these calculations, for completeness, we briefly 
discuss two critical exponents at the transition: v, the correlation length exponent and r/ the anoma- 
lous dimension of the Neel field, n = z*az at criticality. Large- computations for these quanti- 
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ties produce the results^Si^i^, 
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Note that these values become negative and unphysical for the case of interest (A^ = 2), so they do 
not provide any quantitative information for = 2. Presumably, the higher order corrections are 
large and cannot be neglected. It is interesting however to note that the anomalous dimension, r] is 
expected to be large, since at = oo it approaches 1, which is a result of the fact the Neel order 
parameter is not the field that renders the action quadratic at mean-field level. This observation 
agrees qualitatively with the QMC data on the JQ model,^ where a scaling analysis close to the 
phase transition found an anomalous dimension, r] ^ 0.35, that was almost an order of magnitude 
larger than that of the conventional 0(3) universality class. Reassuringly the first term of order 
1/A^ is of the correct sign (negative) correcting the the A^ = oo result, i] = 1, in the correct 
direction. We now turn to large-N computations of certain scaling functions and amplitudes ratios. 



III. FINITE SIZE SCALING FUNCTIONS FOR CF^^ MODEL AT iV = oo 

In this Section we will work only at A^ = oo, but will study arbitrary values of the parameter 
LT. We need to first extremize /° to obtain the large- A^ mass equation in a box: 

,,,2 ,1,2 I 1.2 I ^2 - / ^l2 (^^-^ 



L^Ly ^ ul +kl +kl +ml^^ J 8n^p 

^ nx,ny,nT "-^ 

where Un = 27rTn and kn = and this equation has to be solved self-consistently to obtain the 
saddle point value, mbox- 

We can simplify Eq. (flOl) by using the Poisson summation formula, 

where i = x,y,T and Lr = (3. Taking the N = term on the RHS of Eq. (flOl) . we obtain 
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FIG. 1 : Large- saddle-point value of the mass as a function of the parameter, x = LT. The asymptotes 
for x <C 1 and x S> 1 are plotted for comparison. 



The integral on the LHS can be evaluated exactly in 3 dimensions, giving us the following simple 



non-linear equation for m 



■mho 
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where tn = y {N^ + Ny)x^ + N^, where x = LT. A solution of this non-linear equation is 
shown in Fig. [T] where the self-consistent mass is plotted as a function of the scaling parameter 
X = LT. When x > 1, we know m{x) = 21og(^^^) ^ 0.962424 and for x < 1, we know 
m{x) ^ 1.51196/C, these asymptotes have also been plotted for comparison. 

We now compute the spin stiffness, ps, the uniform susceptibility, Xu, and the specific heat, Cy, 
for the CP^^^ model, when it is placed in a box of linear dimension L and at temperature T. We 
note that due to the absence of any anomalous scaling dimension in ps and Xu, the scaling forms 
are completely universal. In this section we restrict ourselves to the N = oo case. 

In the scaling limit, proximate to the critical point, we can quite generally write. 
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where t measure deviations from the critical coupling, z is the dynamic exponent and c is a non- 
universal velocity. In our large— calculations, we set c = 1 and hence ignore its presence; in our 
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QMC calculations on the other hand, c is determined by the details of our JQ model and we have 
to estimate it from numerical simulations. 

Now that we have calculated the value of the mass parameter, we can calculate the = cx) 
value of Xu and as a function ofx = LT. 
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Completing the Matsubara sum and re- writing in units of T, we find. 
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The sums on and n^, converge fairly rapidly and we can evaluate them by simply introducing 
cutoffs in the sums. We know from simple hyper-scaling laws that Z(x oo, 0) = A^x^. 
In order to evaluate we convert the sums in Eq. (fTSi) into integrals and the RHS becomes 
A^^/N = 0.0856271, the mean-field value of the universal amplitude. 

The computation of ps follows in exactly the same way, though with space and time inter- 
changed. Again we can extract the amplitude, defined by the limiting behavior Y(x — 0, 0) = 
Ap/x, by converting sums into integrals, we find AyN = 0.0926013. The universal functions 
Y(x, 0), Z(a;, 0) can be evaluated numerically and are shown in Fig. |2l 

Comparison with JQ model: The corresponding finite-size scaling functions for ps and Xu 
for the JQ model have been computed before in Ref. 18|. They are reproduced here in Fig. Ob). 
There is a qualitative agreement between the = oo calculations [Fig. [2ta)] and the numerical 
data. It is encouraging that the numerical data also shows the correct asymptotic forms for the 
scaling functions at large LT for Xu and small LT for p^. In order to avoid the complication 
of determining the non-universal velocity, c, it is useful to consider ratios of numbers where this 
quantity cancels. A completely universal number can be constructed by estimating TZ = Ap-J'A^. 
At = oo it takes the value: Tl^=oo = 0.076642. The quoted value of this combination of 
amplitudes in Ref. 8 is 0.075(4). The normalizations of the quantities ps and Xu in the large-N 
and QMC analysis is presented in Appendix lAl The agreement is surprisingly good. We can go 
a step ahead and compare the amplitudes directly. In order to do so, we need an estimate for the 
non-universal velocity. One way to estimate this quantity is to study the data in Fig. [21 Indeed, as 
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FIG. 2: Universal finite-size scaling forms for Xu and pg. (a) for the CP^~^ field theory, evaluated in the 
N = oo limit, as described in the text, (b) from quantum Monte Carlo on the JQ model close to the Neel- 
VBS transition [reproduced from Ref.(8(] . The functions plotted correspond to Y{x, 0) and Z{x, 0). There 
is no scale factor on the y-axis of this plot, but the x-axis has to be scaled properly with a non-universal 
velocity to make a comparison. 



is clear from the study in the Appendix, when LT = c the system is perfectly cubic, and hence 
the two universal functions plotted must be equal, i.e. c is the value of LT when the functions 
cross. By analyzing our data, we find: c = 2.4(3). Using this value of c, the QMC estimates for 
the amplitudes are = 0.23(6) and Ap = 0.15(2) (using the data from Ref. sj), in reasonable 
agreement with the N = oo estimates (A'^ = 0.171 and A^, = 0.185). 

It would be clearly be interesting to verify that the next 1/A^ correction to this quantity is 
actually small. In the next section we study a quantity for which we have succeeded in calculating 
these corrections. 



IV. WILSON RATIO 

In the previous section we restricted ourselves to a = oo calculation. To go beyond a 
qualitative discussion and actually compare numbers, it is essential to include 1/A^ corrections. We 
turn to this computation in the present section. It is also clearly of interest to focus on amplitude 
ratios which do not depend on the non-universal velocity, c. One such ratio is the so-called Wilson 
Ratio, 

W = ^^^ (19) 

Cv Acv 
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where the second equality make use of the scaling forms, Eq. (fT4l) when L ^ 1/T. Note that 
unlike IZ which requires an amplitude from the LT <^ 1 limit, W is a ratio of thermodynamic 
quantities, and is hence accessible even to possible experimental measurements. 

At = oo we can obtain the value the of both amplitudes analytically from the free energy, 
Eq. ([5]), allowing an estimate of the Wilson ratio at = oo: 

Computations of the 1/A^ corrections are rather technical, requiring tedious analytic and nu- 
merical evaluations. The basic calculation has been set up in Sec. Ull The free energy has to be 
computed at next to leading order at finite-T, but in the thermodynamic limit. At order 1/A^, it 
receives contributions from Gaussian fluctuation of both the Lagrange multiplier A and the gauge 
field Afj_. The resulting correction to the free energy, Eq. ([7]), has to be evaluated numerically 
and then numerical derivatives give the specific heat and susceptibility. Explicit details of these 



calculations may be found in Ref. |20|, we will be content with only presenting the results here. 



Al + A]^^ + A'/ AT 0.08562 - 0.02650 + 0.26106 ^^^^^ ^^^^ 
hm Wo(i/N) = -rrr^ ^r. tt^ ~ 1™ { ^ 0.0645 (21) 

N^2 ^^^^^> A'l.^ + A}.\^ + A}f^ N^2 A^(1.8366) - 0.38368 + 2.9928 ^ ^ 

where the superscript indicates A^ = oo values (these contributions are proportional to A^ and 
evaluated at A^ = 2), and the superscript lA and lA indicate the leading 1/A^ correction from the A 
and A^ fields (these contributions have no A^ dependence). It should be noted that the gauge field 
fluctuations do produce rather large corrections to the A^ = oo amplitudes individually (unlike the 
A terms), so it is unclear how reliably they estimate the role of fluctuations. A proper estimate 
likely requires the inclusion of further terms in the expansion. It is re-assuring however that the 
1/A^ corrections do have the correct sign for both A^ and Ap, with respect to the QMC results (see 
Table U). 

Comparison with JQ model: We extract the amplitudes for Cy and Xu by studying the finite 
temperature data on a 128 x 128 system, close to the phase transition in the JQ model, we use 
J/Q = 0.038. Because the specific heat requires a subtraction of two estimators, it turns out to 
be quite noisy, we hence find it preferable to look at the temperature dependence of the average 
energy, which can be measured very accurately. The QMC data and fits to it are shown and 
described in the inset of Fig. |3l From these fits, we can extract the Wilson Ratio. Wqmc = 
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FIG. 3: Quantum critical scaling of the energy, E and the uniform susceptibility, Xu in the thermodynamic 
regime, T l/L. The size has been fixed to L = 128 in these simulations, the data is converged to 
its infinite size limit within our error bars. Fits of E(T) to the form ue + bsT^ and Xu{T) to the from 
+ b^T are shown as dashed lines. The chi-square merit function for these fits per degree of freedom 
(7 DOF in this case) are 1.91 and 1.96 respectively. From these fits we estimate, A^/c^ = 0.041(2) and 
^Cv/c^ = 0.75(5). 

0.055(5). This number is already fairly close to >VAr=oo = 0.0466, the mean-field value. 



V. CONCLUSIONS 

In this paper, we have studied the CP^~^ field theory at large- - the fixed point at = oo is 
known to be stable at large but finite A^. The proof of existence (or non-existense) of a fixed point 
for the case of interest, = 2, is highly non-trivial and does not yet exist. The large- A^ limit, 
however, allows a controlled expansion of universal amplitudes and scaling functions that can be 
extrapolated to the case of interest, A^ = 2. In particular, we have studied a number of amplitudes 
and scaling functions related to finite size and finite temperature effects; two "universal" amplitude 
ratios that can be extracted without the knowledge of c were also studied, TZ = Ap-J'A^ and the 
Wilson Ratio W = A^^j Acv- The Wilson ratio is of greater interest since it is well defined in the 
thermodynamic limit (allowing for instance a possible comparison with experiment). 

As argued in Ref. 11, the low energy behavior close to a transition between the Neel and VBS 
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Definition 


QMC 


large-iV for iV = 2 




limr_^^ y„ = (A^/c^)T 


0.23(6) 


0.17125 


Ap 


limr^o Ps = cApl L 


0.15(2) 


0.18520 


Ac. 


\imL^^C, = {AcJc')T' 


4.3(3) 


3.6733 


W 


AJAc. 


0.55(5) 


0.46622 


n 


Ap \/ Ay^ 


0.075(4) 


0.076642 



TABLE I: Table of amplitudes and amplitude ratios. The first column are the amplitudes (ratios) with error 
bars determined from Quantum Monte-Carlo as detailed in the text. Note that the amplitude ratios, 7^, W 
do not require an estimate of the non-universal velocity, c. The second column is obtained by setting N = 2 
in the = cxd, CP^~^ saddle point theory; these values are known in principle with arbitrary accuracy. 

phases is expected to be described by the non-compact CP^ field theory. It is hence interesting 
to compare our large-N results with the analysis of the JQ model as a test for quantum criticality 
in the JQ model. To facilitate such a comparison, we have estimated a number of the equivalent 
universal numbers mentioned in the previous paragraph from quantum Monte Carlo simulations 
on the JQ model. We have provided a catalogue of our estimates of these amplitudes in Table. HI 
The qualitative agreement for the finite-size scaling functions of Xu and ps in Fig. [2] and the quan- 
titative comparison amplitude ratios studied here is encouraging; all the amplitude ratios at mean 
field agree reasonably with the QMC data. A fully convincing demonstration would require direct 
simulations of an appropriate discretization of the CP^ field theory on large lattices and compar- 
ison with the numerical values we have provided here. This is an exciting direction for future 
work. 
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APPENDIX A: NORMALIZATION OF p, AND Xu 

In this appendix we provide a summary of how we have defined ps and Xu both in the QMC 
calculations and the large-N expansion. The normalization has to be done properly to make a 
numerical comparison between the large-N and quantum Monte Carlo. 

We define Xu as the response to a uniform twist along the temporal direction. 

= ^ (Al) 

where is the angle of the twist along the 2;— direction of the spin per unit of imaginary time and 
fe is the free energy per unit volume calculated with the imposed twist. It is easy to see that this 
definition reproduces the familiar meaning of Xu- 

Vfg = -T In Zg = -T hi Tr [e-^^e^^^"^*'°t ■ ■ ■ ] = -T hi Tr [e-^^e^^'^^''°'] (A2) 
X« = ^ = ^((^tot)^) (A3) 

where we have used the fact that [S^^^., H] = 0. 

Now we can define ps in exactly the same way, as the response to a uniform twist along one 
spatial direction, say x, 

f. - ^ (A4) 

where is the angle of the twist along the 2— direction of the spin per unit length of space and 
is the free energy per unit volume calculated with the imposed twist. It is clear that the twisted 
partition function must be periodic in (pL, i.e. Z^cj)) = J2w Zw^^^'^^ (for exactly the same reason 
that the partition function of charged particles on a ring are periodic in the flux that threads the 
ring), where W is summed on integers and is the winding number of the trajectories of the bosons 
that one would obtain by interpreting our spin model as hard-core bosons. Then by applying the 
formula Eq. (|A4I) . we arrive at the classic result, 

Ps = T{W^), (A5) 

where W is the so-called spatial winding number. For the full original derivation of this idea, see 
Ref. 22! and for an adaptation to the SSE method used here, see Ref. 
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Both these quantities can be calculated by imposing a similar twist to the CP field theory 
which results in modifying either the temporal(spatial) derivative as the case may be, 

drZa drZa + iO-^Zp (A6) 

dxZa dxZa + i(t)-Y^li 

The stiffness or susceptibility is then evaluated as the second derivative of the twisted free energy. 
Such a procedure has been carried out in the body of the text. 
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